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The current-current correlator method has been used successfully to obtain very accurate results 
for quark masses and the coupling a s . The calculations were done using Highly Improved Stag- 
gered Quarks (HISQ) and heavy-heavy meson correlators. We now extend this work to the signif- 
icantly more challenging heavy-light case, reporting the first results here. The aim is to determine 
nonperturbative Z factors for NRQCD heavy-light currents, but first we test the method in the 
HISQ case where Z = 1. 
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1. Motivation 

In the study of semileptonic and leptonic processes, like B — > nlv and B — > Iv, the nonpertur- 
bative Z factors for heavy-light currents are needed. One way to try to calculate these (nonperturba- 
tively on the lattice) is the current-current correlator method, that has been successfully tested and 
used in the heavy onium case [jj], §]. We now want to extend these results and use the same method 
in the heavy-light case. Here we report on the first results using the HISQ action. The eventual aim 
is to extract NRQCD heavy-light Z factors. 



2. Current-current correlator method 

The idea is to match time moments of meson correlators to energy-derivative moments at 
q 2 = of polarization functions n calculated in continuum QCD perturbation theory to high order. 
The pseudoscalar current-current correlators are defined as 

G{t) = a 6 X>m,) 2 (0| 75 (x, f)j 5 (0,0)|0>. (2.1) 

.V 

Then the time moments are 



(see e.g. [jjj ||]). To help reducing the errors we divide each moment by the tree level value, G„\ 
and define reduced moments R n as 



RT = -^y and R^=l-^\ for n > 6. (2.3) 



In the continuum the reduced moments are 

* r = Jp and KT = ^Jy f"r n>6. (2.4, 

The g n are perturbative series in a s (fl), known for the heavy-heavy case through a s 3 (jU) [Ql and for 
heavy-light through ccf(p) [§]. The mass ntf, is the heavy quark mass in the MS scheme at the scale 
fX. Comparing the lattice and continuum R n allows us to extract the mass ratio m^J (2m/ 1 (/i)), and 
thus the quark mass. The calculation above is for the case with no Z factor. If the lattice current 
has a Z factor then that can also be extracted. 



3. Heavy-light JJ correlators 

We compare lattice calculations to continuum perturbation theory through OC 2 (fl). In this 
work we have used coarse, fine, superfine and ultrafine MILC lattice configurations. We calculate 
heavy-light correlators using the HISQ action ^ for both quarks with several heavy quark masses 
from charm up to the b quark mass. Note that Z = 1 in the HISQ case. Some of the calculated 
reduced moments R n are shown in Fig. [j] as examples of our results. Comparing heavy-strange, 
heavy-charm and heavy-heavy correlator reduced moments shows that we can clearly distinguish 
between these three cases. In the following subsections we address some challenges of the heavy- 
light calculations. 
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Figure 1: Top figure: Heavy-strange correlator reduced moments R„ as a function of heavy-heavy meson 
mass Mm (in GeV). The other two figures show the heavy-strange, heavy-charm and heavy-heavy reduced 
moments R4 and R\q as a function of M^. The range is from charm (at about 3 GeV) to b (about 10 GeV). 
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ultrafine heavy-strange, m =0.0165, m =0.5 
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Figure 2: The reduced moments in the free, non-interacting theory depend on the volume. 
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Figure 3: The ratio of reduced moments calculated on two different coarse lattices (one with L = 20, other 
one with L = 28) shows that there is no volume dependence in the interacting case. 



3.1 Volume dependence 

The tree level (free) moments depend on volume — note that this is an artifact of the free 
case only. This is illustrated in Figures || and |3| — the R n depend on volume in the non-interacting 
theory, but not in the interacting case. Therefore we need to calculate the free level moments in the 
infinite volume limit. We do this by calculating the free moments using different volumes, L 3 , and 
fitting them with 

e -A 2 L 

A Q +A X — . (3.1) 
The result in the infinite volume limit is then simply given by the fit parameter Aq. 
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3.2 Quark condensate 

The quark condensate appears in the reduced moment R n at tree level as 
4?r 2 („ _ i) (n _ 2 ) („ - 3) [- a + f ] ( mi W ) 



m 

I | (n-3)mi m 4 
mi, " 



(3.2) 



The quark condensate is not present in the heavy-heavy case, but it is sizeable in the heavy-light 
case — the fraction of tree level qq condensate in R n can easily be 10-30% for heavy quark masses 
masses between c and b, as can be seen in Fig. ^. Note that the leading term is \jrrv\. This 
poses a challenge, as the a s corrections to the condensate are not known. The gluon condensate 
contribution is much smaller and can be safely neglected in the analysis. 

3.3 mi /m-h corrections to perturbative series 

Perturbation theory with m q = is not sufficient, as mi{\l)jmh{\i) corrections become im- 
portant for B c : m c (n)/mt,(iA) « 0.22. At small values of the ratio the m/(/x)/m/j(/i) expansion is 
good, i.e. it works for B s . At large values of the ratio the expansion is not good enough. This is 
illustrated in Fig. ||. However, we now have the exact coefficients for given ratios m/(/x)/m/j(/i) 
for tree level (shown in the plot as bursts) and order a s , so this problem can be partly avoided. The 
exact coefficients are still needed for aj and higher orders. 

3.4 Fits 

We fit the lattice data R]f\ n>6, with 

R n 1 = ( ^, m ?' \ ) ( 1 + c i a s + c 2 af + c 3 q s 3 + c 4 af + c 5 a 5 + C(, af + qq condensate ) 

\2m h (n)J (33) 

(l + b\ (am h (ii)) 2 + b2(amh(n)) 4 + h {am h (\i )) 6 + d\a 2 + d 2 a 4 ) 
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HL moment coeffs - exact vs expansion in ml/mh HL moment coeffs - exact vs expansion in ml/mh 




Figure 5: m/(/i)/m/,(jU) corrections to continuum perturbation theory coefficients. Here x — mi/m/,. The 
tree level values of the coefficients have been divided by the tree level value at mi /mi, = 0. 



and choose the scale /i = nth- We take the first few coefficients (ci and in the heavy-light case) 
from continuum perturbation theory, and treat the coefficients for the higher order a s terms as fit 
parameters. The quark condensate is given in Eq. 32 at tree level. We take the qq condensate 
value to be (m s ss) = (0.2 GeV) 4 from the Gell-Mann - Oakes - Renner relation, allowing the ss 
condensate to be 0.7 times the light quark condensate. We also allow for the presence of higher 
order condensate terms estimating them with powers of the leading condensate. In B c there is 
no condensate contribution, and we get a good fit using the exact coefficients (order a s ). As the 
lattice calculations are done at a non-zero lattice spacing a, we include a-dependent terms in the fit 
function — even powers of a and am^(^i). 

To extract the mass ratio m^ A /(2m/ 1 (/i)) we use the lattice simulation data (Eq. 2.3), with 
am rjA / anih from the lattice simulations, and compare these /?Jf tt to the continuum perturbation the- 
ory result (Eq. 2.4). That is, we find values for o^jg(ju) and /n rj;i /(2m/,(/x)) that make lattice and 
continuum results agree for small n > 4. This can then be combined with experimental results for 
the rib, v\c meson masses to obtain the quark masses. In the heavy-light case we can use the oiy^{n) 
values extracted from the heavy-heavy calculation. 

To test the method in the heavy-light case we look at the mass ratio m, jA /(2m/ ! ( i u)) and com- 
pare to heavy-heavy results. This is shown in Figure ^. In the heavy-strange case the fits are to 
one R„ at a time, not to all R n simultaneously as in the heavy-heavy case. The mass ratio extracted 
from the heavy-strange correlator moments is the same as in the heavy-heavy case, as expected, 
but currently a lot less accurate. 



4. Conclusions and future 

We are extending the use of current-current correlator method, earlier used to study heavy- 
heavy systems, to heavy-light systems. The full analysis of heavy-light data is still in progress, 
but we can already say that the JJ correlator method works well. The quark condensate sets some 
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23456789 10 
Figure 6: The mass ratio m IJ/i /(2m/ I (/i)). 



limitations — we can not use high moments in the B s case. However, in B c there is no condensate 
contribution. Our aim is to extract Z for NRQCD — there the challenge is to control relativistic 
corrections. 



References 

[1] C. McNeile, C.T.H. Davies, E. Follana, K. Hornbostel, G.P. Lepage; HPQCD collaboration; 
Phys.Rev.D82:0345 12,2010, arXiv: 1004 .4285 

[2] I. Allison, E. Dalgic, C.T.H. Davies, E. Follana, R.R. Horgan, K. Hornbostel, G.P. Lepage, C. 

McNeile, J. Shigemitsu, H. Trottier, R.M. Woloshyn, K.G. Chetyrkin, J.H. Kuhn, M. Steinhauser, C. 
Sturm; HPQCD collaboration; Phys.Rev.D78:0545 13,2008, arXiv: 0805 .2999 

[3] K. G Chetyrkin, J. H. Kuhn and C. Sturm; Eur.Phys.J. C48 (2006) 107, arXiv : hep-ph/0 604234 

[4] K.G. Chetyrkin, M. Steinhauser; Eur.Phys.J. C21 (2001) 319-338, arXiv : hep-ph/0 10 8 017 

[5] E. Follana, Q. Mason, C. Davies, K. Hornbostel, G.P. Lepage, J. Shigemitsu, H. Trottier, and K. 
Wong; HPQCD collaboration; Phys.Rev.D75:054502,2007, arXiv : hep-lat /0 6100 92 

[6] L. J. Reinders, H. R. Rubinstein and S. Yazaki; Phys. Lett. B 97 (1980) 257-260 



7 



